Origin of the insulating state in honeycomb iridates and rhodates 
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A burning question in the emerging field of spin-orbit driven insulating iridates, such as Na2lr03 
and Li2lr03 is whether the observed insulating state should be classified as a Mott-Hubbard insulator 
derived from a half-filled relativistic j^b = 1/2 band or as a band insulator where the gap is assisted 
by spin-orbit interaction, or Coulomb correlations, or both. The difference between these two 
interpretations is that only for the former, strong spin-orbit coupling (A > W, where W is the band 
width) is essential. We have synthesized the isostructural and isoelectronic Li2Rh03 and report its 
electrical resistivity and magnetic susceptibility. Remarkably it shows insulating behavior together 
with fiuctuating effective S = 1/2 moments, similar to Na2lr03 and Li2fr03, although in Rh*"'' 
(4d"') the spin-orbit coupling is greatly reduced. We show that this behavior has non-relativistic 
one-electron origin (although Coulomb correlations assist in opening the gap), and can be traced 
down to formation of quasi- molecular orbitals, similar to those in Na2lr03. 

PACS numbers: 75.10.-b,75.10.Jm,71.70.Ej,71.15.Mb 



In recent years complex iridium oxides have caused ex- 
traordinary interestplSl since the physics there is governed 
by a unique combination of several comparable scales: 
one-electron hopping t, spin-orbit coupling (SOC) A, and 
the Hubbard repulsion U. The honeycomb layered com- 
pound Na2lr03, a small-gap antiferromagnetic insula- 
^QjEEl^ has received particular attention. It was suggested 
that the adequate description of the electronic behav- 
ior of this system is in terms of a half-filled relativistic 
Jeff = 1/2 band, which becomes a Mott-Hubbard insula- 
toi^l. However, U in iridates is rather small (1.5-2 eV), 
therefore the corresponding band must be rather narrow 
for the system to become insulating. In this scenario, 
the SOC is the leading interaction in these systems, so 
that the t2g bands split into a narrow doublet with the 
effective angular moment jcS = 1/2 and a quartet with 
Jeff = 3/2. The one-electron hopping is suppressed by a 
factor of 1.5A/t ^ 3, where 1.5A is the energy separation 
between the doublet and the quartet. 

Recently, some of us^ proposed an alternative descrip- 
tion and argued that the one-electron non-relativistic 
band structure might be a better starting point for the 
description of the electronic behavior of honeycomb iri- 
dates than the limit X ^ W (band width). In this case 
the band structure is dominated by the formation of so- 
called quasi-molecular orbitals (QMOs), and consists of 
four narrow bands (the width being defined by second- 
neighbor hoppings and other secondary one-electron pa- 
rameters), spread over a width of ~ At, where t ^ 0.3 
eV is the leading one-electron hopping. The highest and 
the lowest bands are singlets, having one state per Ir (i.e., 
one state per spin per unit cell of two formula units), and 
the two middle bands are doublets. In the A = limit the 
upper singlet and doublet bands nearly merge, forming 
a triplet manifold, while turning on SOC further splits 
those bands into three singlets. The upper two bands 
barely overlap, forming an incipient (SOC assisted) band 



insulator, and even a very small C/ of a few tenth of an cV 
is sufficient to open a gap. In this picture, the material is 
characterized as a spin-orbit assisted insulator with the 
gap enhanced by Hubbard correlation. 

In view of the two alternative, and partly oppo- 
site, descriptions of the insulating state in these sys- 
tems, we present here a comparative analysis between 
the electronic behavior of hexagonal iridates and rho- 
dates. Specifically, we have synthesized and investigated 
Li2Rh03, which shows insulating behavior at low tem- 
peratures, similar to Na2lr03 and Li2lr03, even though 
in Rh'^+ (M^) the SOC is substantially reduced. The 
comparison sheds light onto the nature of the insulating 
state in these systems. 

The paper is organized as follows. We first settle the 
terminology between the various definitions of insulators. 
We then proceed with electrical resistivity and magnetic 
susceptibility data of Li2Rh03 and the description of its 
electronic and magnetic properties by means of density 
functional theory (DFT) calculations with and without 
inclusion of spin-orbit coupling and discuss the similari- 
ties and differences of this rhodate system compared to 
the hexagonal iridates. As an outlook, we provide some 
predictions for the magnetism in the hexagonal rhodates. 

The question of whether a particular phase is charac- 
terized as a Mott-Hubbard insulator or a band insulator 
is largely terminological, as no strict definition or crite- 
rion exists to rigorously separate these notions. Some 
authorgSl further subdivide the classification of insulators 
into Peierls, Wilson, Slater or Hund insulators, to men- 
tion a few. We feel that this fine tuning is not helpful 
here, and will concentrate on the difference between Mott 
and band insulators, which is fundamental in the sense 
that one cannot go from the former to the latter continu- 
ously. Note that this division docs not have a one-to-one 
correspondence with the strongly correlated - weakly cor- 
related dichotomy; a band insulator, in our terminology. 



may have a gap strongly enhanced by correlations, but is 
still "topologically connected" with an uncorrelated in- 
sulator. 

We can illustrate this with a simple example: imag- 
ine a crystal of atoms with one half-occupied orbital. 
If the crystal has one atom per unit cell, then on the 
one-electron level this material can never be insulating. 
Upon including an onsite Hubbard U, penalizing double 
occupati on an d hindering itinerancy, it becomes a Mott 
insulator^^^, with no coherent quasiparticles. This hap- 
pens roughly when U becomes larger than the total band 
width W. Now, suppose the same atoms are bound in 
dimers forming a molecular crystal. Each dimer devel- 
ops a bonding and an antibonding state split by some 
energy A. We now allow inter-dimer hopping. The levels 
will broaden into bonding and antibonding bands of the 
width W. As long asW<A, the material is a band insu- 
lator. If A ~ VF, the gap is very small, and the indirect 
gap may even become negative. If we add a Hubbard U to 
this system (not necessarily larger than W) the gap will 
get larger by some fraction of U (depending on the de- 
gree of itinerancy) and this may be a substantial enhance- 
ment. We call this a correlation- enhanced band insulator. 
For instance, solid Ne is a band insulator, even though 
in LDA calculations i ts g ap is severely underestimated 
(12.7 eV vs. 21.4 eV).'^^ This discrepancy is related to 
another problem in the density functional theory (DFT) , 
the so-called density derivative discontinuity, and not to 
Hubbard correlations. To first approximation, this dis- 
crepancy is inversely proportional to the static dielectric 
function.ESl 

An example of a Mott insulator is FeO. It has one 
electron in the spin-minority t2g band, and is a metal in 
DFT. Coulomb correlations have to destroy entirely the 
coherent DFT metallic band crossing the Fermi level, and 
the excitation gap appears between the incoherent lower 
and upper Hubbard band.I^ Note that despite FeO be- 
ing a Mott insulator even in the paramagnetic phase, the 
simplistic treatment of LDA-|-U,^ as opposed to more 
sophisticated DMFT,^^ cannot reproduce insulating be- 
havior by symmetry; the cubic symmetry needs to be 
broken, for instance by assuming antiferromagnetic or- 
dering, after which a gap opens at sufficiently large val- 
ues of U. Similarly, in parent compounds of the super- 
conducting cuprates there exists one hole in the Cg band, 
and symmetry does not allow to open a gap in DFT. 
These systems are Mott-Hubbard insulators. 

Finally, MnO is an example of a (strongly) correlation- 
enhanced band insulator. It has a gap between 3d major- 
ity and minority bands. In DFT, this gap is driven by the 
Hund's rule, and is ^ 51 — W, where the Stoner factor / is 
~ 1 eV. The calculated value is 1.4 eV as compared to 4.5 
eV in the experiment.— This material is strongly affected 
by Mott physics, and routinely called Mott insulator, yet 
one can make a gedanken experiment and gradually re- 
duce the Hubbard correlations to zero, whereupon the 
gap will drop to its DFT value, without losing the in- 
sulating character. Note also that a Mott-Hubbard in- 



sulator, in our nomenclature, does not necessarily imply 
a strong Hubbard repulsion U ^ t, where i is a typical 
intersite hopping. For instance, TaS2 by no means can 
be expected to be a strongly correlated material, and U 
cannot be more than a fraction of an eV, and, indeed, at 
high temperatures it is a metal. Yet at low temperature 
it experiences a charge density wave transition typical 
for this structure, which, combined with the spin-orbit 
interaction on Ta, splits off, essentially accidentally, an 
ultra-narrow band {W '--^ 0.1 eV), and even a minuscule 
U suffices to split it into two Hubbard bands. 

In Na2lr03, U is relatively small and the material can- 
not be strongly correlated. Moreover, correlations are 
suppressed in LDA-I-U calculations because of substan- 
tial delocalization of electrons over Irg hexagons, so that 
in order to increase the calculated gap from ^ to ~ 0.3 
eV one has to add U ^ A eV. Furthermore, there is in- 
direct evidence of itinerancy in the experiment: the or- 
dered magnetic moment even at the lowest temperature is 
less than 0.3/iB,l3in reasonable agreement with the band 
structure calculations (w 0.5 /is, equally distributed be- 
tween spin and orbital moments), while the fully local- 
ized jeff model (l/z_B, split 2:1 between spin and orbital 
moments) requires strong fluctuations to suppress the or- 
dered moment .121 

It is often argued that the experimentally measurecP^ 
spin-orbit correlation factor, (L • S) , is consistent with 
one hole in the jeff = 1/2 state, and thus proves its ex- 
istence. However, this factor is mostly collected from 
the Cg holes ,t'^-'^-^ and is well described by band structure 
calculations. 

A comparison of honeycomb iridates with the isostruc- 
tural and isoelectronic Li2Rh03 should be very instruc- 
tive, because if the former are SOC Mott insulators like 
Sr2lr04,El^then a Rh analogue (with much weaker SOC) 
should be metallic, just as Sr2Rh04.^- If, on the other 
hand, the formation of quasi-molecular orbitals triggers 
the insulating behavior, then a larger U in Li2Rh03 will 
likely recreate the same physics as for Na2lr03, i.e. a 
correlation-enhanced band insulator. 

We have synthesized Li2Rh03 polycrystals by the solid 
state reaction method from stoichiometric amounts of 
Li2C03 and Rh powder. The mixture has several times 
been pelletized and reacted in O2 flow at temperatures up 
to 850°C. Powder X-ray diffraction (XRD) scans do not 
reveal any evidence for secondary phases and are simi- 
lar to those reported in Ref. [24j see Fig. [2] Magnetic 
susceptibility and (four-probe) electrical resistivity have 
been determined utilizing commercial (Quantum Design) 
instruments. 

As presented in Fig. [I] (a), Li2Rh03 shows clear in- 
sulating resistivity behavior, which follows the same 
variable-ra nge h opping dependence as found in Na2lr03 
or Li2lr03.'^j25] Previous resistivity measurements at 
higher temperatures found an activation gap of ^ 
80 meV.f^S The magnetic susceptibility (Fig. [l] (b)) is 
Curie- Weiss (CW) like, with a small kink at 5 K, likely 
due to some spin-glass freezing, which needs to be inves- 




TABLE I: Comparison of the first principles hopping amph- 
tudes in Li2Rh03 and Na2lr03. The notations are explained 
in detail in Ref. i2lj 



100 200 300 50 100 150 200 250 300 
7(K) 



FIG. 1: Temperature dependence of the electrical resistiv- 
ity (a) and magnetic susceptibihty (b) of polycrystalhne 
Li2Rh03. Inset (c) displays the resistivity data (on log 
scale) vs. T"!/*. The hne in (b) is a Curie- Weiss (CW) fit 
X{T) = Xo + C/{T - Qw with xo = -1.235 x 10"* cm^mol 
and Qw ~ —59 K. The inset (d) shows Ax~^ vs T with 
Ax = xiT) - Xo, and the CW fit. 



tigated in future measurements. The CW fit between 100 
and 300 K corresponds to ii^s — 2.2/.tB. Similar results 
have been recently reported by Luo et al:-- 

Now, we need to establish the crystal structure. Lab 
powder XRD is insensitive to the 0-positions, which hin- 
ders structural determination. For instance, initial re- 
finement for Na2lr03 suggested substantial Na/Ir site 
exchange,'^ but later measurements on a single crystal, 
showed that the leading defects were stacking faults of 
well ordered honeycomb planes.^^' Similar problems are 
expected for Li2Rh03.^'* For that reason, we have used 
(well-determined) unit-cell parameters for Li2Rh03, and 
performed first principles optimization of the internal 
parameters.^ We note the same procedure yielded ex- 
cellent agreement with the refined crystal structure of 
Na2lr03. The final structure is presented in Table \n\ 
This refined structure is consistent with the lab powder 
XRD data (Fig. [2]). 

Despite the overall low crystal symmetry, the lo- 
cal symmetry of the Rh2Li planes is rather high: the 
hexagons are nearly ideal and the Rh-O-Rh angles are 
nearly the same and relatively close to 90°. This makes 
it a showcase for the quasi-molecular orbital concept.'^ 
To this end, we have performed first principle calcula- 
tions using the WIEN2k code,^^ and projected the results 
using a standard Wannier function projection technique 
as proposed by Aichhorn et alP^ and further developed 
in Ref. 1291 The resulting tight-binding parameters are 
shown in Table H] 

As we see, the main condition for the QMO picture 
(dominance of the 0-assisted nearest neighbor hoppings) 
is fulfilled. Projecting the density of states (DOS) onto 
individual QMOs we see that, although it does not sepa- 
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*For the shorter distance, the first number corresponds to 
xy — > xz, xy —^ yz transitions and the second to xz —> yz 

transitions. 

^For the shorter distance, the first number corresponds to 

xy —> xy transitions and the second number to xz —> xz, yz —> yz 

transitions. 



rate into isolated manifolds as in Na2lr03, it is composed 
of overlapping QMOs as shown in Fig. |3] 

We have also performed spin-polarized calculations 
with various spin configurations'^^ (Fig. H(a)). We were 
not able to stabilize a Neel order (magnetic moments 
collapse), but the ferromagnetic (FM) and two antifer- 
romagnetic phases, the "stripy" and the "zigzag" phase, 
are all stable, with the ground state practically degener- 
ate between the two AFM states. The FM state has a 
small advantage in the calculations, which is lost upon 
application of U (see below). One has to keep in mind 
though that DFT has the tendency to overstabilize ferro- 
magnetism in oxides, due to overestimation of the Hund's 
coupling on oxygen. The calculated FM state, just as in 
Na2lr03, is a half-metal with M = l/i^/Rh. 

Neither the nonmagnetic state (Fig. Isl) nor any of the 
magnetic states considered (FM, stripy and zigzag) are 
insulating. Including SOC has little effect on either en- 
ergetics or proximity to an insulator (Figs. HIs] and l6|. 

On the other hand, experimentally this material ap- 
pears to be insulating. It is natural to attribute this fact 
to the effect of Hubbard U, which in 4d metals is about 
3-4 eV, twice as large as for 5d systems. Even though 
LDA-|-(on-site)U is a rather naive way to tackle correla- 
tions in a QMO system, we have tried, faut de mieux, 
to a pply a standard LDA-fU correction to our calcula- 
tions.'^ As expected, for C/ > 3 eV we obtain an insu- 




FIG. 2: (Color online) Comparison of observed (black) and calculated (red) powder XRD spectra for Li2Rli03 (see text). 




FIG. 3: Non-relativistic nonmagnetic band structure and 
density of states of Li2Rh03, projected onto quasi-molecular 
orbitals, as described in Ref. 1211 



I 



> -0.1 

CD 



'-0.2 
-0.3 
-0.4 



Q[\ 



> -0.1 

CD 



-0.2 



-0.3 



▲ ;;;;7^ 

nonmagnetic ■■■" 

stripe ■■♦■■ 

zigzag ■♦■ 

ferromagnetic ■■^■■ 

(a) without SO coupling 



0.4 - (b) with SO coupling 



2 

U-J (eV) 



lator for both AFM configurations. In Figs. [5] and |6] we 
show the evolution of the density of states (DOS) with 
the Hubbard U for the zigzag and stripy configurations, 
respectively. Besides, adding U produces somewhat less 
obvious effects. First, it destabilizes the FM structure, 
making all three magnetic structures degenerate within 
the computational accuracy (in calculations with SOC 
the FM state is a few meV lower in energy, but, as men- 
tioned, DFT always slightly overestimates this tendency 



FIG. 4: Energy of different magnetic configurations in eV/Rh 
relative to the nonmagnetic state, as a function of (U — J). 
Energies at ([/ - J) / are offset by 1.05{U - J). The top 
panel (a) does not include spin-orbit coupling, the bottom 
panel (b) does. 



because it includes Hund's rule self- interaction on O). 
Second, the Hubbard U enhances the SOC, increasing 
the calculated orbital moments. The spin moment also 
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TABLE II: Optimized crystal structure of Li2Rh03, using 
experimental lattice parameters (a = 5.123 A, 6 = 8.836 A, 
c = 5.885 A, /3 = 125.374°) and space group C2/m. 
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FIG. 5: Evolution of the density of states with the Hubbard U 
in the zigzag phase without spin orbit coupling (left panels) 
and with spin orbit coupling (right panels). We use Jh = 
0.7 eV throughout. 
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FIG. 6: Evolution of the density of states with the Hubbard U 
in the stripy phase, without spin orbit coupling (left panels), 
and with spin orbit coupling (right panels). We use Jh = 
0.7 eV throughout. 



positively correlates with U, but the dependence is much 
weaker. For instance, for the FM state the moment inside 
the Rh muffin tin sphere increases from 0.58 to 0.66/iB 
as U increases from to 5 eV. 

An important point to make is that, as one can ex- 



pect from the small value of the SOC, it is not essential 
for obtaining an insulating state; an antiferromagnctic 
order, however, is, just as in such prototype Mott insu- 
lators as FeO and CoO. In the recent literature one can 
find an opposite statement, that SOC is a requirement for 
reproducing the insulating behavior.'^Sl This stems from 
a misconception that LDA-I-U is capable of describing 
Mott insulators in the paramagnetic case. One of the 
pathologies that LDA+U shares with LDA is absence of 
local magnetic fluctuations. In both methods instead of 
a paramagnetic state, i. e. a state with disordered local 
moments, a non-magnetic state, with no moments at all 
is considered. As a result, in such prototype strongly 
correlated materials as for instance 3d oxides a metal- 
lic state is protected by symmetry, unless some magnetic 
ordering is included (in some cases even the ferromag- 
netic order suffices, in others an antiferromagnctic order- 
ing is needed), and LDA-I-U fails to reproduce the para- 
magnetic insulating phase. Li2Rh03 is a similar case. 
In the nonmagnetic calculations there are band cross- 
ings protected by symmetry, that is to say, the best one 
can possibly achieve within LDA+U, even with an arbi- 
trarily large U, is a zero-gap semiconductor. SOC, even 
infinitesimally small, removes this protection (the pro- 
tected bands can now hybridize), and now a sufficiently 
large U can open a full gap. Obviously, this fact does not 
tell us anything about the real role of the SOC, but only 
highlights shortcomings of the LDA+U method.'^ 

The fact that three different magnetic configurations, 
FM, zigzag and stripy, with ordered moments on Rh 
(~0.5-0.7 /xs) independent of the magnetic pattern, are 
very close in energy indicates considerable frustration. 
Structural disorder is then expected to push the system 
towards a spin-glass regime. 

These results show an important similarity between 
the 5d compound Na2lr03 and the isostructural and iso- 
electronic Ad compound Li2Rh03, despite a much larger 
Hubbard U and much smaller spin-orbit A in the latter. 



This similarity suggests that properties of these materials 
are largely controlled by the non-relativistic, one-electron 
physics, namely the formation of quasi-molecular or- 
bitals, while the role of Coulomb correlations and SOC 
lies primarily in enhancing already existing tendencies 
(in particular, toward insulating behavior). As a word 
of caution, we want to emphasize that while our results 
point toward these systems being band (Slater) insula- 
tors rather than Mott insulators, this does not indicate 
that they are weakly correlated or that they are local- 
ized rather than itinerant. On the other hand, our results 
suggest that local antiferromagnetism is an important in- 



gredient in the formation of an insulating state and that 
Coulomb correlations are instrumental in enhancing the 
insulating gap. 
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